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The relativistic generalization of the Newtonian Lagrangian perturbation theory is investigated. 
In previous works, the first-order trace solutions that are generated by the spatially projected gravi¬ 
toelectric part of the Weyl tensor were given together with extensions and applications for accessing 
the nonperturbative regime. We furnish here construction rules to obtain from Newtonian solutions 
the gravitoelectric class of relativistic solutions, for which we give the complete perturbation and 
solution schemes at any order of the perturbations. By construction, these schemes generalize the 
complete hierarchy of solutions of the Newtonian Lagrangian perturbation theory. 

PACS numbers: 98.80.-k, 98.80.Jk, 04.20-q, 04.20.Cv, 04.25.Nx 


I. INTRODUCTION 

In previous work of this series of papers we laid down 
the foundations of the Lagrangian perturbation theory 
by writing Einstein’s equations in 3 + 1 form for a sin¬ 
gle dynamical variable. We investigated its first-order 
solutions for the trace and antisymmetric parts, we ex¬ 
trapolated this solution in the spirit of Zel’dovich’s ap¬ 
proximation in Newtonian cosmology, and we provided a 
definition of a nonperturbative scheme of structure for¬ 
mation [LI]. We then studied the average properties of 
the latter in relation to the Dark Energy and Dark Mat¬ 
ter problems in Ref. [L2] . Here, we proceed by providing 
the gravitoelectric subclass of relativistic nth-order per¬ 
turbation and solution schemes. As in previous work we 
restrict our attention to irrotational dust continua for 
simplicity. The generalization to more general matter 
models is scheduled. 

The problem of perturbation solutions in general rela¬ 
tivity (GR) has been addressed by a plethora of works. 
In cosmology the ‘standard approach’ is based on the 
gauge invariant ‘Bardeen formalism’ (for a selection of 
key references on standard perturbation theory see Refs. 
[4], [35], [28], [19]). A covariant and gauge-invariant ap¬ 
proach has been proposed [22, 23], together with var¬ 
ious other approaches not listed here (some being dis¬ 
cussed within a variational framework in a recent paper 
[30]). The reason for the existence of various approaches 
is due to an ambiguity of the choice of perturbation vari¬ 
able, the choice of a ‘background’, but also due to dif¬ 
ferent philosophies; e.g. the standard gauge-invariant 
approach compares the physical manifold with a refer¬ 
ence ‘background manifold’, while others solely operate 
on the physical manifold. The conceptual difference of 


our framework lies in the fact that we no longer consider 
a reference background manifold. All the quantities are 
now defined on the physical space section. All orders 
of the perturbations are defined on the physical mani¬ 
fold, not with respect to a zero-order manifold (that was 
interpreted as the background manifold in standard per¬ 
turbation theory). Moreover, we are perturbing a single 
dynamical variable which, intuitively, is the square root 
of the spatial metric using the ‘Cartan formalism’. As 
a consequence, the issue of gauge invariance does not 
arise; covariance or diffeomorphism invariance is guaran¬ 
teed for a given foliation of spacetime by using Cartan 
differential forms. We shall address the representation 
of the perturbations in other foliations of spacetime in a 
followup article. 

A similar point of view has also been taken in previous 
work, i.e. the pioneering work by Kasai presents a rel¬ 
ativistic generalization of the ‘Zel’dovich approximation’ 
[46], and followup works with his collaborators present 
a class of second-order perturbation solutions [27, 41]; 
see also the earlier papers by Tomita [43-45], the paper 
by Salopek et al. [42] as well as the series of papers by 
Matarrese, Pantano and Saez [32-34], considerations of 
so-called ‘silent universe models’ [5, 6, 8, 24, 29], and the 
recent paper [38]. These works are all in a wider sense 
concerned with the relativistic Lagrangian perturbation 
theory and concentrate on an intrinsic, covariant descrip¬ 
tion of perturbations. Still, the present work takes an¬ 
other angle and goes beyond some concepts of these latter 
works through the following elements: 

• We consider, as in Refs. [LI, L2], a formalism that 
allows us to write the Einstein equations within 
a flow-orthogonal foliation with a single dynami¬ 
cal variable comprising the spatial Cartan coframe 
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fields. These furnish the conceptual generalization 
of the Lagrangian deformation gradient being the 
single dynamical variable in the Newtonian theory. 
One advantage of this approach is that only per¬ 
turbations of this variable are considered, which 
entitles us to express all other physical quantities 
as functionals of this variable. Thus, it is possible 
to leave the strictly perturbative framework and to 
construct nonperturbative models by injecting the 
deformation solutions at a given order of expansion 
of the Einstein equations into the functional defini¬ 
tions of these fields, without a posteriori expanding 
the functional expressions. This in turn provides 
highly nonlinear approximations for structure for¬ 
mation (e.g., the density field is known through an 
exact integral of the perturbation variable; the met¬ 
ric as a bilinear form maintains its role as a measure 
of distance, i.e. as a quadratic expression; the cur¬ 
vatures are the general defining functionals for the 
given perturbed space, etc.); 

• We provide construction rules to derive relativis¬ 
tic perturbative solutions from the known Newto¬ 
nian solutions at any order of the perturbations: we 
have to additionally study the traceless symmetric 
part of the equations having no obvious Newtonian 
analog, and which is fundamentally linked to the 
traceless Ricci tensor and the physics of gravita¬ 
tional waves. In the present work, however, we re¬ 
strict our attention to that subclass of the traceless 
symmetric perturbations that are generated by the 
spatially projected gravitoelectric part of the Weyl 
tensor. We show that this part, in turn, can be 
constructed from the tidal tensor of the Newtonian 
theory. For this purpose we employ a division of the 
governing equations into gravitoelectric and grav- 
itomagnetic parts; 

• We give the perturbation and solution schemes to 
any order of the perturbations for the gravitoelec¬ 
tric part of the Lagrange-Einstein system. These 
schemes cover the full Newtonian hierarchy of the 
Lagrangian perturbation theory using a restriction 
rule that we will define. This allows us to construct 
the leading-order modes of relativistic solutions at 
any order. Thus, the results of this work allow 
us to construct the relativistic counterpart of all 
Newtonian results, where higher-order information 
is needed, e.g., to construct the bispectrum of the 
perturbations (cf. Refs. [37, 40]). 

• We perform a strictly intrinsic derivation, i.e., with¬ 
out reference to an external background space. The 
perturbations are described locally (in coordinates 
of the tangent spaces at the physical manifold at 
each order of the perturbations). There is no need 
for a diffeomorphism to a global background man¬ 
ifold that also, in general, does not exist; existence 
of a global diffeomorphism depends on the global 


topology of the physical manifold which, in gen¬ 
eral, needs several coordinate charts to cover it. 
Note that in the literature on Lagrangian relativis¬ 
tic perturbations mentioned above, although start¬ 
ing with the Cartan formalism, the nonintegrabil- 
ity of the Cartan deformations is given up for the 
building of solutions, hence implicitly introducing 
a reference background space for the perturbations. 

Before we start, let us recall our strategy (for details 
the reader is directed to Ref. [LI]). In the Newtonian 
theory the Lagrangian picture of fluid motion allows us 
to represent Newton’s equations in terms of a single dy¬ 
namical variable, the Lagrangian deformation gradient 
built from the trajectory field. For this system the gen¬ 
eral perturbation and solution schemes at any order are 
provided in Ref. [20], a paper to which we henceforth 
often refer. 

Einstein’s equations within a flow-orthogonal foliation 
of spacetime can be formulated in terms of equations for 
the gravitoelectric and gravitomagnetic parts of the spa¬ 
tially projected Weyl tensor. Subjecting the gravitoelec¬ 
tric subsystem of equations to a “Minkowski Restriction ”, 
i.e. by sending the Cartan coframes to exact forms, we 
obtain the Newtonian system in Lagrangian form [LI], 
[14]:Sect. 7.1. In this paper we investigate the reverse 
process, i.e. the transposition from integrable to non- 
integrable deformations, which enables us to construct 
a gravitoelectric subclass of the relativistic perturbation 
and solution schemes that corresponds to the Newtonian 
perturbation and solution schemes. 

While the Newtonian system furnishes a vector the¬ 
ory, where the gravitational field strength is determined 
by its divergence and its curl (the trace and antisymmet¬ 
ric parts of the Eulerian field strength gradient), the so 
generalized schemes deliver nontrivial solutions for the 
trace-free symmetric part that is connected to the gravi¬ 
toelectric part of the spatially projected Weyl tensor, the 
Newtonian counterpart of which is the tidal field tensor. 

The paper is structured as follows. Section II recalls 
the equations of Newtonian and relativistic cosmology 
for an irrotational dust matter model. We highlight a 
formal correspondence between the Newtonian equations 
and the relativistic gravitoelectric part of the equations 
by employing a geometrical restriction procedure, named 
Minkowski Restriction. In Section III we investigate per¬ 
turbation and solution schemes at any order n of the 
perturbations by explicitly paraphrasing the Newtonian 
schemes. Section IV explains the reconstruction rules 
and provides explicit examples. Finally, Section V sums 
up and discusses perspectives. 


II. EQUATIONS OF MOTION AND 
CONSTRAINTS 

In this section, after setting notations, we recall the 
Einstein equations, written in 3 -I- 1 form and expressed 
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through a single dynamical variable, represented by Car- 
tan coframe fields as functions of local coordinates in the 
3—hypersurfaces. This recalls the parts of Ref. [LI] rel¬ 
evant to this paper. 


A. Notations and Technicalities 


We employ the differential forms formalism for its com¬ 
pactness and antisymmetric properties and its spatial dif- 
feomorphism invariance. We also project to the common 
coefficient formalism in which we work out the solutions. 

We consider a set of a differential k —forms k“. The 
coefficients of these fields can be expressed in the ex¬ 
act basis {dX 1 } of the cotangent space at a given point, 
k a = ik dX n A ... A dX lk 7 where A is the wedge 
product, the antisymmetrization of the tensorial prod¬ 
uct A/\B = A®B — B® A. Their exterior derivative 
yields dk° = k\ ik<[ip dX h A ... A dX ik A dXV For 
general forms we choose the letters a,b,c-- ■ as counter 
indices (they refer to the nonexact basis), while the let¬ 
ters are reserved for coordinate indices (they 

refer to the exact basis). The Hodge dual is denoted by 
a star and defined in TV— dimensional space by: 


*k“ 


Vy ^ h—ik ,h..-4 
(N — fc)! ik+i—is 


dP^ 1 A ... A dX jN 


with j N the Levi-Civita pseudotensor. 

In most perturbation approaches the bilinear metric 
form is considered as the dynamical variable. In this 
article, we consider the matter model “irrotational dust” 
and employ a 3+1 flow-orthogonal foliation of spacetime, 
for which the 4- and 3-metric bilinear forms read: 


*- 4 ^g = — dt ® d t+ with * 3) g = gij dJ 1 <g> dX- 7 , (1) 


where X 1 are Gaussian normal coordinates, here equiv¬ 
alent to the Newtonian Lagrangian coordinates. The 
resulting split of the system of equations (Arnowitt- 
Deser-Misner (ADM) system) is composed of 6 equa¬ 
tions of motion and 4 constraint equations. In this fo¬ 
liation the 4 Cartan one-forms can be restricted to a 
t— parametrization of 3 spatial one-form fields. 

In general relativity a spatial description of the fluid 
continuum in terms of vector-valued trajectories is im¬ 
possible, unless we move to a higher-dimensional em¬ 
bedding vector space. To describe the fluid intrinsically 
(i.e., without reference to an embedding vector space), 
it is necessary to introduce nonexact forms, known as 
the Cartan spatial coframe fields r/® = g a i dX' 1 , with 
a = 1 • • • 3. The Cartan formalism permits switching 
between a nonexact basis and the coordinate basis. A 
key element is the freedom of choice of the normaliza¬ 
tion of the nonexact basis. In order to obtain equations 
that are formally closer to the Newtonian ones, we do not 
choose orthonormal (Cartan) coframes f) a as is common 
in the literature, but more general ones r/ a that we will 
call adapted coframes. The reader is directed to Refs. 


[LI, L2] for additional information and implications re¬ 
lated to this choice. Formally, this means that the spatial 
metric form is decomposed as: 

(3) g = G a bT] a ® rj b , (2) 


where G a b is constant in time: G a b = G 0 b(X). Note that, 
if coframes become exact forms 77“ = d /“, the counter 
indices become coordinate indices, since the functions /“ 
can be used to define global coordinates x 1 = / a_H . In 
this case the metric can be brought (by a spatial diffeo- 
morphism) into the form: 


(3) g = Sijdx 1 ® dor 7 , (3) 


which defines a flat spacetime (see the proof in ap¬ 
pendix A). 

The exact functional for the density is given as in the 
Newtonian approach: gj = pi, where the index i marks 
the initial conditions and J is defined as coefficient func¬ 
tion of the 3—volume form, normalized by the determi¬ 
nant of the initial metric: 



( 4 ) 


with y/g d 3 X the 3—volume form on the exact basis, 
g := det(gij(X.,t)) and G := det (Gij) = det(f/y(X, ti)). 
We have the relation: 

A r 1 b A r) c = ^JdX 1 A dW A dX fc , (5) 

6 6 


where J = det ( 77 “). 

Below, we first recall the basic systems of equations 
governing an irrotational dust continuum in the La¬ 
grangian formulation of the Newtonian theory. Then, 
after presenting Einstein’s theory formulated in the La¬ 
grangian frame, we list the counterpart of the gravito- 
electric subsystem of equations in the latter theory. 


B. Newtonian Theory 

In the Lagrangian picture of self -gravitating fluids a 
family of trajectories, x l = f l (X k ,t), labelled by their 
Lagrangian coordinates { X l }, i = 1,2,3, is introduced 
[13,17]. It furnishes a one-parameter family of diffeomor- 
phisms, parametrized by the Newtonian time t, between 
the Eulerian, {a; 1 }, and the Lagrangian coordinates. Reg¬ 
ular solutions of the Lagrange-Newton system of equa¬ 
tions have to obey four evolution equations. The three 
components of the trajectory field (also position field) 
f(X k ,t) are the only dynamical variables. Other fields 
are conceived to be represented as functionals of the tra¬ 
jectory field like the velocity and acceleration fields, the 
density and vorticity fields, etc., 

«*:= f, a 1 :=/*, 0 = <?,( J/ J;)" 1 , 


( 6 ) 
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where the overdot denotes time-derivative along the tra¬ 
jectories. J = det(/*| fc ) is the Jacobian of the coordinate 
transformation and J\ is the initial Jacobian, with spatial 
derivatives with respect to Lagrangian coordinates being 
abbreviated by a vertical slash |. The acceleration field 
a 1 is identified with the gravitational field strength g l , 
respecting the equivalence of inertial and gravitational 
mass. Once a given field is represented as a functional of 
the deformation field, it can be written in the Eulerian 
frame by inserting the inverse of the transformation f l . 
Note that Ji can be set to 1 if we require x l = X 1 at 
initial time. The closed Lagrange-Newton system is de¬ 
fined by the nonlinear gravitational evolution equations 
(7) and (8) for the deformation gradient, see Ref. [20]: 

Sijdf* A cl/' = 0 , (7) 

\tijkdP A d.f A d/ fc = (A - 4tt Gg) d 3 / , (8) 

with A the cosmological constant, G the gravitational 
constant, Eq. (7) corresponding to the three field equa¬ 
tions for the Eulerian curl, and Eq. (8) corresponding 
to the field equation for the Eulerian divergence of the 
gravitational field strength. In the above equations the 
exact integral for the density (third equation of Eq. (6)) 
has to be inserted to reduce the number of variables. For 
A = 0 the system does not explicitly contain the Jaco¬ 
bian, provided J / 0. Regular solutions are characterized 
by J > 0. Since, in the Newtonian theory, the Cartan 
coframe fields are exact forms, Eq. (4) reads: 

J : = A df j A df k = ^JdX i A dX j A dX fc 

6 6 

<=> J := Jd 3 X = d 3 /; J = det(/%) , (9) 

where d 3 X is the Lagrangian volume 3—form. The coef¬ 
ficients of these equations are equivalent to the following 
coefficient equations for the deformation gradient d/ ? in 
Lagrangian coordinates: 

Ww = 0 - (10) 

\e ijk e^r b f\ q f\ = A J - 4 ttGqJ . (11) 

An alternative to express Eq. (10) reads: 

S kr e pq[ p ilrn f r] \Af q \ m = °- ( 12 ) 

Expressed in terms of the Newtonian tidal tensor, 

= ^e o6c e ifei /y W |l - J ( A - 51 j ; (!3) 

%]=0 ; S\ = 0, ( 14 ) 

Eqs. (14) correspond to the Lagrange-Newton system of 
equations {(10), (11)}. 

C. Einstein’s Equations in Lagrangian form 

We formulate here the Einstein equations in terms of 
Cartan coframe fields as they are transported along the 


flow lines (here spacetime geodesics). The system of 
equations we obtain will be called the Lagrange-Einstein 
system. 

In terms of coframe fields, the irrotational dust contin¬ 
uum is governed by the following evolution and constraint 


equations: 

Gab rf A rf = 0 ; (15) 

7j £ dbc (ri a A T] b A T] c )' 

= (-H a d + (4nG Q + A)5 a d )Jd 3 X ] (16) 

e abc rf A rf A r) c = (16t tGq + 2A - 77) Jd 3 AT ; (17) 
e abc (dfj a A rf + u) a d Arf A rj b ^j = 0 , (18) 


where the equations are, respectively, the irrotationality 
condition on the gravitational field (15), the equation of 
motion (16), the energy constraint (17) and the momen¬ 
tum constraints (18). For irrotational matter flows, as is 
assumed throughout this paper, the first equation can be 
replaced by the kinematical irrotationality condition: 

Gab r/ a A r/ b = 0 . (19) 

Nevertheless, we consider the double time-derivative ex¬ 
pression for two reasons: (i) in the Newtonian limit we 
want to reproduce the field equations, which involve a 
second time-derivative, and (ii) for a general system, this 
equation is always true because of the conservation of the 
vorticity 2—form, us = G a b rf A rf = oq; for the Newto¬ 
nian case, u) = d (rqdA' 1 ) = uq (see Appendix B in Ref. 
[20]). Nevertheless, it is important to note that the 3+1 
foliation formalism cannot describe a nonzero vorticity. 
The description of a nonzero vorticity will require a 1 + 3 
threading of spacetime (see, e.g., Ref. [25]) . 

The combination of the trace of the equation of motion 
and the energy constraint straightforwardly leads to the 
Raychaudhuri equation: 

7 fabcV a A rf A rf = (A - AttGq) Jd 3 X . (20) 

To derive the above equations we have implicitly used the 
Cartan connection one-form and the curvature two-form 
that we do not need explicitly in what follows: 

u a b ■= l a cbV c , ( 21 ) 

n a b :=^n a bcdV c Arj d , (22) 

with the connection and curvature coefficients r y a cb and 
77 a bcd in the nonexact basis, respectively. The 3-Ricci 
tensor can be expressed through the curvature two-form: 

n a d rf Arf Arf = 5 db VL a d A rj° - S dc n a d A rf . (23) 
We employ the Hodge star operation to obtain the coef- 
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ficient equations in the exact basis dX l \ 

Gab V[i r l b j] = 0 ; (24) 

jjtabce' kl (fi^vWi)' = -Ki + ^nGg + A) 5^ ; (25) 

^e«6ce mj V TO >7 6 i»7 c fc = ~ + (8*Gg + A) ; (26) 

{eabce^vWi)^ = > ( 27 ) 

where a double vertical slash denotes the covariant spa¬ 
tial derivative with respect to the 3-metric and the spa¬ 
tial connection is assumed symmetric. As before, the first 
equation can be replaced by the irrotationality condition: 

GabV[iV b j] = 0. (28) 

Again, the trace of the equation of motion and the energy 
constraint leads to the Raychaudhuri equation: 

^je abc e m 7i a iV b k V c e = A - 4 nGg . (29) 

The system {(24) — (27)} consists of 13 equations, 
where the first corresponds to the irrotationality condi¬ 
tion (3 equations), the second to the symmetric evolution 
equations (6 equations), subjected to 4 constraint equa¬ 
tions that are the ADM constraints (1 equation for the 
energy constraint, and 3 equations for the momentum 
constraints). Thus, the first 9 equations furnish evolution 
equations for the 9 coefficient functions of the 3 Cartan 
coframe fields. 

The above system is equivalent to the results devel¬ 
oped in Ref. [LI] in a different basis: in the first paper 
the choice of the standard orthonormal coframes has been 
made, whereas since Ref. [L2] the choice of the adapted 
coframes is preferred for reasons of allowing us to con¬ 
struct a formally closer Newtonian analogy. 

D. Equivalence of the two gravitoelectric sets of 
equations in the Minkowski Restriction 

We will now discuss the link between parts of the rela¬ 
tivistic system and the full Newtonian system. Formally, 
this link is provided by the Minkowski Restriction. 

1. Definition of the Minkowski Restriction 

Let T) a be Cartan one-form fields in a 4—dimensional 
manifold (Greek letters are used in 4 dimensions). A 
set of forms r] a is said to be exact , if there exist func¬ 
tions /“ such that r) a = d/“, where d denotes the exte¬ 
rior derivative operator, acting on forms and functions. 
The Minkowski Restriction (henceforth MR) consists of 
the replacement of the nonintegrable coefficients by in- 
tegrable ones, r] a v —r f keeping the speed of light 
c finite. With this restriction, the Cartan coframe coef¬ 
ficients yield the Newtonian deformation gradient, and 


the local tangent spaces all become identical and form 
the global Minkowski spacetime. The Newtonian limit 
could be defined as the MR of Einstein’s theory and ad¬ 
ditionally sending c to infinity. In the flow-orthogonal 
foliation, employed in this paper, the 4—dimensional 
coframes reduce to r/ a = (dt,r/ a ), and their MR reads 
d/“ = (df,d/ a_H ). Note that c and the signature are 
carried by the 4—dimensional metric coefficients; c is set 
to 1 throughout this paper. We will need the inverse MR 
and use it as a rule to construct relativistic Lagrangian 
solutions from known Newtonian solutions. (For the MR 
of the metric, compare appendix A, and the remark on 
the Newton-GR “dictionary” in appendix B). 


2. Gravitoelectric equations 

Reference [LI] noted that a part of Einstein’s equa¬ 
tions, namely {(15), (20)}, are related to the gravitoelec¬ 
tric part of the spatially projected Weyl tensor. This 
tensor is trace-free and represented by the 3 one-form 
fields E“ (see Ref. [LI] Eq. (A23)): 

E a = -r, a + l(A-4nGg)r ] a -, (30) 

then, the irrotationality condition (15) and the trace 
equation of motion (20) are generated by 

G ab E“Ar, 6 = 0 ; e abc E Q A V b A rf = 0 . (31) 

These two equations are therefore referred to as the grav¬ 
itoelectric part of Einstein’s equations. A projection of 
the gravitoelectric one-form fields and Equations (31), 
using the Hodge star operator, yields to their coefficient 
representation: 

= -^jeau^rjVWi + l ( A - 4 ttGq) 6^ ; (32) 
E m = 0 ; E\ = 0. (33) 

(Note that Eij := 8 b i GbaE a 'j = G a iE “.) This just pro¬ 
vides a rewriting of {(24), (29)}. (A remark on the grav- 
itomagnetic part of the spatially projected Weyl tensor 
can be found in appendix B.) 


3. Executing the MR 

Sending the spatial Cartan coframes to exact forms, 
i.e., executing the MR, their coefficients rf' i are restricted 
to the Newtonian deformation gradient f °L The gravito¬ 
electric system of equations in the form of Eqs. (33) then 
reduces to the Newtonian system in the form of Eqs. (14); 
note the conventional sign change between the gravito¬ 
electric part of the spatially projected Weyl tensor E l - 
and the Newtonian tidal tensor £® •. This operation closes 
the system, reducing the number of free functions from 
nine \q a i {X k ,t)) to three (f 1 (X k ,t)). A consideration of 
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the MR for the remaining equations, yielding nontrivial 
Newtonian analogs, will not be needed in this paper, but 
will be the subject of forthcoming work. 

Considering only the gravitoelectric equations is not 
enough to determine the nine functions of the coframe 
coefficients. The relativistic aspects contained in the re¬ 
maining gravitomagnetic equations will lead to a richer 
structure of the solutions and also to constraints on so¬ 
lutions of the gravitoelectric system. A followup work 
will explicitly consider both parts in the framework of 
first-order solutions. 

To conclude: the Lagrange-Einstein gravitoelectric 
equations are (up to nonintegrability) equivalent to their 
Newtonian analogs, whereas the gravitomagnetic equa¬ 
tions have no obvious Newtonian counterpart (this issue 
is more subtle than the mere absence of a gravitomag¬ 
netic analogy, as we explain in the followup paper). 


III. CONSTRUCTION SCHEMES FOR 
RELATIVISTIC PERTURBATIONS AND 
SOLUTIONS AT ANY ORDER 


We now turn to the main part of this paper and con¬ 
struct the gravitoelectric subclass of nth-order relativis¬ 
tic perturbation and solution schemes through general¬ 
ization of the known Newtonian schemes. This allows 
furnishing relativistic inhomogeneous models for large- 
scale structure formation in the Universe. The successful 
Lagrangian perturbation theory in Newtonian cosmology 
is well-developed. We will generalize here the perturba¬ 
tion and solution schemes of Newtonian cosmology given 
in the review [20] , the essential steps of which will be re¬ 
called in this section, followed by their relativistic coun¬ 
terparts. 

All schemes are applied to the matter model ‘irrota- 
tional dust’. It is possible to extend the present schemes 
by employing the framework for more general fluids in a 
Lagrangian description that will be developed in forth¬ 
coming work. Most of the known representations are 
focused on writing equations in terms of tensor or form 
coefficients. Our investigation will be guided by the com¬ 
pact differential forms formalism as before. However, we 
will also project to the coefficient form in parallel to ease 
reading. 


A. General nth order perturbation scheme 


As in standard perturbation theories, we decompose 
the perturbed quantity into a Friedmann-Lemaitre- 
Robertson-Walker (FLRW) solution and deviations 
thereof, which are expanded up to a chosen order n of 
the perturbations. Contrary to the standard perturba¬ 
tion theory, we do not perturb the metric globally at the 
background space, but we perturb the Cartan coframes 


locally: 


rf = rf i dX l = a{t) ( d° + ]T dA* , (34) 


in the local exact basis dA*. Notice that with this ansatz 
we choose to perturb a zero-curvature FLRW model, but 
it is possible to encode an initial first-order constant cur¬ 
vature in the coefficient functions G a b in the following lo¬ 
cal metric coefficients, which can be calculated from the 
above coframe ansatz: 


9ij = GabV a iV b j ■ (35) 

Furthermore, we can link these results to the ones ob¬ 
tained for the orthonormal coframes f) c (compare also 
the corresponding remarks in [18]: p.35 and [26]: p.75). 
Indeed, the metric bilinear form can be written as: 


g = S cd f) c ®f] d = G ab V a ®V b ■ (36) 


From this identity, we conclude: 

Gab = S cd rfj] d b , (37) 

where the fj c a are the coefficients of the projection of 
rf onto the basis rf. In the next subsection, we will 
specify the coframes we consider in such a way that the 
initial coframe perturbations vanish. From now on, we 
will call these coframes adapted coframes to distinguish 
them from the orthonormal ones (see also appendix A). 


B. Initial data for the perturbation scheme 

We choose initial data in formal correspondence with 
the Lagrangian theory in Newtonian cosmology and gen¬ 
eralize these initial fields to the relativistic stage. This 
has obvious advantages with regard to the aim to give 
construction rules that translate the known Newtonian 
solutions to general relativity. For the initial data set¬ 
ting in the Newtonian case, see Ref. [20]. 


1. Fundamental initial data 

Let the three one-form fields U° = [/“dA 1 be the ini¬ 
tial one-form generalization of the Newtonian peculiar 
velocity-gradient, obtained by the inverse MR. Accord¬ 
ingly, let W“ = W“ dA* be the initial one-form general¬ 
ization of the Newtonian peculiar-acceleration gradient. 
Our solutions will be written in terms of these initial 
data. They determine the initial values of the coframes 
as follows: 

V a i(t i) = <5“ , (38) 

(ti) = Hid" + U a i ; Hi := cq ; 0l := 1 ; (39) 

ri a i (t i )=a i 6 a i + W" i . ( 40 ) 
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Equation (38) implies that the coframes we will work 
with from now on are initially equal to the exact La- 
grangian coordinate basis: 77 “ (tj) = dX l . This in 
turn provides the initial metric coefficients in the form: 

G ij =g ij (t i ) = G ah 5 a i 5 b j . (41) 

In view of the flow-orthogonal foliation, we have the 
irrotationality constraint: 

ui = G ab r] a A rj b = 0 =► G ab U a A S b j dX j = 0 . (42) 

This implies for the coefficient functions: = 0. (We 

used the implicit definition Uij := S^G^U ^.) 

Remark: 

From (37) and (38), it is interesting to notice the follow¬ 
ing relations that hold to zeroth and first order (the full 
initial data are considered to be first order, as was also 
the choice in the Newtonian schemes [20]): 

r G® = Sir, 

J G$= 2 #y; (43) 

i 2 Ab)=G«+2P fe) ; 

where = P,; 7 - (t{). We are thus able to rederive some 
results from the ones obtained in previous works that 
used orthonormal coframes. For example, the Ricci cur¬ 
vature tensor at first-order can be obtained by injecting 
the identities (43) into (93) of Ref. [LI]. We can so 
obtain the adapted coframes from the orthonormal ones 
and vice-versa (compare also appendix A). 


2. Relativistic counterpart of the Poisson equation and 
consequences for W“ 

In the Newtonian approach the initial peculiar- 
acceleration and the density inhomogeneities are linked 
through the Poisson equation. To generalize this equa¬ 
tion to the relativistic case, we note the relativistic gen¬ 
eralization of the Newtonian field strength gradient that 
follows from inspection of the Lagrange-Einstein system 
(for details the reader can always consult Ref. [LI]): 

r 3 := 0V + e^e*- 

= -TVj - 00 ' ? + ( 4 tt Gq - A) 6 i j + Q i k Q k j , ( 44 ) 


In terms of the coframe fields, the relativistic gravita¬ 
tional field can be written as follows: 

Fj = ^je abc e M v a jV \v c i ■ (46) 

(For an alternative derivation using a Newton-GR “dic¬ 
tionary” see appendix B). 

Hence, inserting the coframe perturbations and evalu¬ 
ating this expression at initial time, we get the following 
relations (note that the zero-order fields trivially satisfy 
the second constraint): 

T k k (U) = A - AnG Qi = 

A — 47rGgjji(l + 5i) = 3hi + S k a W a k ; 

= S^GbaW'x = w [ij} = 0 . , (47) 

with the initial density contrast <5i. Thus, the deviation 
one-form fields W a obey the following equations that 
generalize the Poisson equation for the inhomogeneous 
deviations off the zero-order solution: 

* U abc W a A S\dXi A S c k dX k = -AirGSg-, ; 

2 J 

G ab W a A 6 b j dX j = 0 , (48) 

with 5qi = Qi — QHi: implying for the coefficient functions: 

-^ k a W a k = 8 ei = em 5 i ; W M = 0. (49) 

3. Summary of initial data 

We summarize the set of initial data, determined by 
our choice of the basis and subjected to the constraints. 
We assume in perturbative expansions, without loss of 
generality [20], that the initial data {(50) — (52)} are 
first order. We drop the index (1) for notational ease 
and denote the initial data for the comoving perturbation 
form coefficients by P“(tj) =: &. We set: 

• for the initial deformation and the initial general¬ 
izations of the Newtonian velocity and acceleration 
gradients: 

= 0 \j n ^ 

U a(1) = U“, U {ij] = 0 ; (50) 

W a(D = W a ; ^ = 0 ^ 


with the 3—Ricci tensor coefficients IZij of which the 
trace is the Ricci scalar 1Z , and 0jj the expansion ten¬ 
sor coefficients. According to the energy constraint, 
TZ + 0 2 — <d k e Q e k = 6nGg + 2A, the symmetry of the ex¬ 
pansion tensor and Ricci curvature, it is straightforward 
to show that the relativistic gravitational field coefficients 
Fij respect the following field equations: 

F\ = A - 4t xGg ; F m = 0 . (45) 


• where the coefficients are related via the initial val¬ 
ues of the time-derivatives of the deformation: 


= U a i ; 

& a . = w a i - 2H i U a i 


(51) 


• together with additional initial constraints that are 
to be respected (a relation to the initial metric, 



to the initial density contrast, and the four ADM 
constraint equations evaluated at initial time): 


Gij = Gab^it OTJ^i) = GabS^Sh 
AttGSq^ = -W ; 
nit,) 


H-JJ = -- 


— w 


(52) 




C. Gravitoelectric Perturbation Scheme 

We now recall the general Lagrangian perturbation 
scheme of Newtonian cosmology and generalize it to a 
gravitoelectric scheme in relativistic cosmology. By con¬ 
struction, this latter will already contain the known La¬ 
grangian perturbation scheme at any order in the geo¬ 
metrical limit of exact deformation one-forms. 


(here and in the following we use the abbreviations 
S k a U a k =: U, S k a W a k =: W for the trace expres¬ 
sions.) 

• the initial Ricci curvature as found from the equa¬ 
tion of motion (25): 

n){U) = -(Wj + HiUj) - 5 i j (W + HiU) 

—£ ab c£ ilk U a jU b i5 c k ; (53) 


1. Recap: Newtonian Theory 

The general perturbation scheme has been fully devel¬ 
oped in Ref. [20]. Our approach only slightly differs in 
terms of the initial conditions: we formulate them such 
that they are formally closer to the relativistic approach. 
Following the general ansatz (34), we introduce three co¬ 
moving perturbation forms dP 1 of the three components 
of the comoving vector perturbation fields P l {X l ,t): 


• equating this expression with the initial Ricci ten¬ 
sor as calculated from the initial metric, 


■n ' \ _ 1) /~ib\a(l) . b|a(l) i(l) 

jRiJ — ( - T [j| 61 tjr a +Lr i Wt’k 


I /^a(i)/^[i>|i](i) I 1 /~Ia,i) b\i(i) _ 1 ) 

+ U b\j + 2 Lr {a\b} U j 2 Gq ^b'l 6 ] 

1 

a *~i[b\i](l) O/^C 1 ) i/^[b\ a ]W 

/(jr b ^ [j | o.l ( - T Mbl > 


_|_ 1QW a Q<- 1 ') b \i 


2qW kl4 
pk] 


2G (2) 

p'k] 


T m b] 


(54) 


where, e.g., G ia = G vj 8' J a and G l 3 = G kj S k \ we 
determine the first-order part of the initial metric 
(which is a derived quantity), 


2G ( ^J fcl = -H,UT - W) - ( H,U + W) 5U , (55) 

as well as the second-order part of the initial met¬ 
ric (which later appears in the perturbation and 
solution schemes): 


2G 


(2) [i|fc] 

im 


W^w)) 


(56) 


where the function / can again be derived by equating 
(54) and (53). All further order right-hand sides of (53) 
vanish. 

The initial data given in (50) are exhaustive: in our 
ADM split, the system of equations {(24) — (27)} con¬ 
tains 9 second-order equations of motion for the coframes 
subjected to 4 constraints. A general solution therefore 
contains 18 coefficient functions encoded in Uij and W,j 
that reduce to 12 functions for solutions of the irrotation- 
ality conditions (24), the latter of which are represented 
by the 6 constraints U[ij ] = 0 and = 0. The general 
solution is further subjected to the 4 ADM constraints 
resulting in corresponding constraints on Uij and W,j. 


df(X,t) =■■ a(t)dF*(X,t) = a(t) (dX* + dP\X, t)) , 

(57) 

and decompose the perturbation gradient field on the 
FLRW background order by order: 


dP i = Y e m d P i(m) . (58) 

m= 1 


It is, of course, possible to consider perturbations of the 
position fields f l , because the Newtonian equation can 
be expressed in a vectorial form. The relativistic equa¬ 
tions are, however, tensorial, and we therefore consider 
the representation in terms of the gradient of the fluid’s 
deformation. 

To provide unique solutions of the Newtonian system, 
suitable boundary conditions have to be imposed. For 
the cosmological framework the requirement of periodic 
boundary conditions for field deviations from a Hubble 
flow is a possible choice [16]. This translates into an 
integral constraint on the perturbations: integration over 
a compact spatial domain A4 implies the following: 


[ dP* = f P i = 0 

M dM 


oo 

pi _ f rapi(m) 

m —1 


(59) 


Recall now that IT = dW = U) 3 -dX j and W* = dW = 
W l \jdXi are the initial one-form peculiar-velocity gra¬ 
dient and the initial one-form peculiar-acceleration gra¬ 
dient. The fields W l are determined nonlocally by the 
following set of equations, equivalent to Poisson’s equa¬ 
tion: 


W\i = *^e ijk dW z A dX j A dX fc = -4 tt G6qi ; 

SijdW* A dX 1 = *d (rdl 4 ) = 0 . (60) 
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In view of the restriction to irrotational flows, we addi¬ 
tionally impose the constraint: 

Sijdf A df 3 = 0 => SijdU* A dX j = *d (f/’dl 1 ) = 0 . 

(61) 

Without loss of generality, we can choose the following 
general set of initial data that can be obtained in the 
Newtonian theory or else from the MR of ((50)-(52)): 

• for the initial deformation, peculiar-velocity and 
peculiar-acceleration: 

d^ l(n) = 0 Vn ; 

dU i(1) = dlP, U m = 0 ; (62) 

dW iw = dr, W m = 0 , 

• together with the definition of the Lagrangian met¬ 
ric coefficients and the initial data relation to the 
density perturbation: 

c k c i 

1 (63) 


9ij — $klf |if |j ; 
dpi (1) = Sg, = gmSi = - 


TT^ wk \k ■ 

47 tG ]k 


The metric is Euclidean, since the coefficients can be 
transformed to the coefficients Sij with the help of the 
to / inverse coordinate transformation. 

Plugging the ansatz (57) into the Newtonian equations 
{(7), (8)}, we find for the background Friedmar 
tion: 


e iik 3—dX i A dX j A dX k = 
a 

e ijk (A - AttGqh) dX l A dX j A dX k 

==> 3— = A — 47 tGqh , 
a 

and a full hierarchy of the perturbation equations: 
Sij dP i A (dX 3 + d P j ) = Sijd - 2 dip A dX> ; 


tijk 


(2?id P l ) A dP A dX k + (2V 2 dP i ) A dP 3 A dX k 


+ (V 3 d P l ) A dP j A dP fc 
4 t tG 


— Gjk - 


Sgia^dX 1 A dP A dX k 


where we defined the operator 


De '-^ + 2H Jt ~ T Gqh 


( 66 ) 

(67) 


(The reader may note a difference in the numerical coef¬ 
ficients from the Ref. [20], which we had to correct, see 
appendix B.) Projecting with the Hodge star operator to 
the coefficient form (and integrating Eq. (65)), we obtain: 

p m = JPmWiP^dt' ; (68) 

(69) 


T>iP\i = -AnGSgia 
1 




Imn 


p\, p 3 V 3 P k 


1 1 r \r. 


'3-r \ n 1- 2d | 


P 3 t VoP k 


\n 


After splitting Eqs. (65) and (66) order by order, we 
obtain n sets of equations. At first-order we get: 

(70) 

(71) 

(72) 


SijdP im A dX j = 0 ; 
f-ijkDi dP i(1) A dX j A dX fc = 
a- 3 e ijk dW i A dX j A dX fe ; 
in coefficient form: 


p (1 ) _ n 

Mdi] ~ U 


DTi™ = a~ 3 WV 


'U |i -« •'*' | i 

i.e., a set of linear equations. The generic nth-order sys¬ 
tem of equations will be written below with an implicit 
summation over the order of perturbations in the source 
terms: 

A (P) B <9) = Y A w B {q) , (73) 

p+q—n 

A M B M C m = Y A (r) B^C (t) . (74) 

r-\-s-\-t=n 

Thus, at any order n > 1, the perturbation equations 
read: 

SijdP A dX 3 = -SijdP iw A dP jM ; (75) 
eij k T>idP i( - n) A dX j A dX fc = 


^ijk 


(2V 2 dP Hp) ) A dP 3M A dX* 


equa- 

+ (V 3 d P i(r) ) A dP jM 

A dP fc(t) 

1 

(76) 


in coefficient form: 





p(") _ f p(p) pm(q) Jj 

m ~ Jt m|[ili] 



(77) 

(64) 

v piw = _l imnpjWptWp pi(r) 

1 \i 2 \ rn \ n ° r 


(65) 

- [v 2 P % |f) P 3 \f 

+ (v 2 P 

(p)A 
\j ) 

Pf. (78) 


The reader may consult the review [20] and references 
therein for further details. 


2. Einstein’s Theory 

Assuming the perturbation ansatz (34) for the 
coframes, and using the operator T>i as defined in (67), 
the analogous expansion is performed: the zeroth order 
again leads to the Friedmann equation, and the general 
perturbation scheme reads: 


G ab P a A 6 b jdX 3 + G ab P a A P - 0 ; 
e a bc ViP a a S b jdX 3 A 5\dX k 

+ {2V 2 V a ) A P b A S\dX k + ( V 3 P a ) A P b A P c 

= e abc wl ; a- 3 S a i dX i A S b q dX 3 A S c k dX k . 

3 J 


(79) 


(80) 













10 


In coefficient form and integrating Eq. (79) they become: 
P [tj] = G ab £ ifiPhfldt'; (81) 

V X P\ = - ((2? 2 P\) P 3 j - (V 2 P)) P 3 ) 
~e ljk e lmn (p 3 P ! |; ) P\ m P\ n + Wa~ 3 . (82) 


The hierarchy begins with the first-order equations 
{(70), (71)} which are uniquely determined by the con¬ 
straint initial data (62). The general nth-order, n > 1, 
solution scheme from Eqs. {(65), (66)}, reads: 

SijdP i(n) A dX 3 = N S {n) ; (90) 

djkVidP^ A dX 3 A dX k = N T (n) , (91) 


Expansion order by order leads to the first-order gravi- 
toelectric equations: 


and the general nth-order, n > 1, set of nonlinear equa¬ 
tions: 


uniquely determined by the source terms: 


G ( °l P° (1) A dX b = 0 ; 

(83) 

jl 

e abc PiP a(1) A S b jdX 3 A 8\dX k 
r 3 e ab c W Q A S b j dX i A 5\dX k ; 

(84) 

N q~(n) _ 


f 

Jto 


dpi(P) A d pi^dt' 


(92) 


+ (P 3 dP l(r) ) A d P 3(b) A dP fc(f) l .(93) 


G ab P Q<n) A S b jdX j = -G a b P aM A P b(9) ; (85) 

eabcP iP“ (n) A 8 b j dX 3 A 5%dX k = 

-e abc [2 (V 2 P a(p) ) A P b(9) A 5 c k dX k (86) 
+ (I»3P aW )AP'' w AP c( ‘ ) ] . 

In coefficient form, this reads: 

P$ = 0 ; = Wa" 3 , (87) 

and 

P M= G % £ P V Wpb i\W ■, (88) 

ViPT = -\e ijk e lmn P 3 £P k ™ (p 3 P} m ) 

- {v 2 P i( i P)S ) P j - q) + (P 2 P*j p) ) P j i q) ■ (89) 

This provides equations for the perturbation fields at any 
order n from solutions of order n — 1. 

Comparing {(85), (86)} to the Newtonian equations 
{(75), (76)}, we see (not surprisingly) that we arrive 
at two equivalent sets of equations if we link the per¬ 
turbations P“ and dP* at any order via the MR - 
recall that the construction was done by inversion of 
the MR. dP i = P i ]j . dX j P“. dX 3 = P“; for the 

initial data: d IT = U\ dX 3 ^ [/“ dX j = U a and 
d W l = W l \jdX 3 i—>■ W a jdX 3 = W a . Therefore, we can 
simply translate the formal solution scheme for the trace- 
parts and the antisymmetric parts of the perturbations. 
However, note already here that the inversion of the MR 
produces a symmetric traceless component that is repre¬ 
sented in Newtonian theory by the tidal tensor. 


We have earlier demonstrated the formal equivalence be¬ 
tween the Newtonian equations and the relativistic grav- 
itoelectric equations. The generalization of the Newto¬ 
nian solution scheme to obtain the corresponding rela¬ 
tivistic scheme is now straightforward. 


2. Einstein’s Theory 

The perturbative gravitoelectric Lagrange-Einstein 
system starts at n = 1 with the Eqs. {(83), (84)}, 
uniquely determined by the corresponding constraint ini¬ 
tial data (50). The nth-order, n > 1, gravitoelectric 
solution scheme reads: 

G ab P a(n) A 5 b jdX 3 = <S (n) ; (94) 

e afcc PiP a(n) A 8 b j dX i A 5%dX k = T (n) , (95) 

which is uniquely determined by the source terms: 

S M := G% f (-P a(s) A P b(t) ) dt' ; (96) 

T M := -e a6c (2 (V 2 P aM ) A P bM A 5 c k dX k 

+ (P 3 P aM ) A P b(3) A P c(t) ) . (97) 

The coefficient form of these equations is given by Eqs. 
(87)-(89). 


IV. APPLICATION OF THE SOLUTION 
SCHEME 


D. Gravitoelectric Solution Scheme 

1. Recap: Newtonian Theory 

We first recall the general solution scheme given in Ref. 
[20] , written for the perturbation gradients only. 


To illustrate the use of the scheme {(94), (95)} in prac¬ 
tice, we will in what follows explicitly explain the con¬ 
struction of relativistic solutions from Newtonian ones 
for the general procedure and through examples, in Sub¬ 
sections IV C and IV D, respectively. Before we do so we 
explain the general systematics of the solution scheme. 
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A. Systematics of the solutions 

The nth-order scheme is a hierarchy of ordinary 
second-order differential equations, sourced by an in¬ 
homogeneity resulting from combinations of lower-order 
terms. Thanks to the linearity of the ordinary differen¬ 
tial equations (ODEs), the solution is, at any order n, a 
linear superposition of modes that we will label by l: 

piw = pi.(*’') . (98) 

l 

In the Newtonian case, and for the gravitoelectric rela¬ 
tivistic part, the modes can be further separated into spa¬ 
tial and temporal parts: PL (n,li = ^ n,l \t)P l ^ n ’ l) (Xk). 
This is due to the fact that (95) is an ODE and that its 
coefficients only depend on time. 

From the theory of second-order ODEs it is known 
(see, e.g., Section 2.1.1 of Ref. [36]), that an equation of 
the form 


h (°) y" + /i («) y' + fo (a) y = g (a) , (99) 


will have as the general solution: 

a 

y (a) = Cij/i (a) + C 2 y 2 (a) + J G (a, s) g (s) ds , (100) 

a i 


where Green’s function G (a, s) is defined by 
1 V 2 (a) Vi (s) - yi (a) y 2 (s) 


G (a, s) = 


h (s) 2/i (s) 2/2 (s) - y[ (s) y 2 (s) ' 


( 101 ) 


Therefore, at any order, the solution will have two modes 
l that are given by the homogeneous solution, known for 
a given background model (in the examples we will ex¬ 
plicitly give the solutions for the Einstein-de Sitter case, 
henceforth EdS, and the Cold Dark Matter background 
with a cosmological constant, henceforth ACDM). The 
different modes of the particular solution can be calcu¬ 
lated from the integral in Eq. (100) by setting g = T (n) . 
As integration is linear, the particular solution can be 
computed for each subpart of the source separately, and 
those parts appear as a PL (n,i) in the sum (106). 

To study these subparts, we split the perturbations 
into their trace, their symmetric trace-free part and their 
antisymmetric part: 


pa = pS a dX j 

3 J 


n Q + qr. 

Then, Eqs. {(88), (89)} read: 

” [\ ( G l[i n 1] )ib<S) - G a[Alf P(i) ) 
+ l( G a[^1] pM - G a[^7 P ^ 


W = 


( /^( 7- )or> a ( s )TT^( t ) /^( r )m b ( t )TT°( s ) 

' y^ab V [i 11 j} ^ab [< 11 j} 


) 


( 102 ) 


+ 


; (103) 


Vl pM = _lpMp> 2 pM + U b< a q) V 2 U a l P) 
o 


tp7 ) p 2 f p a 


[pop' 1 -) - ir^n 6 ^ - t pT ) ^ f 'a r 


+ - (-p(‘)<p a w - q3“WpW) p 3 qj^) 

3 

+ - - n a < () p (r >) p 3 n^ s) 

3 

+ (n Q } t) n c ' , ' ) + n c ^n^ r) + n a ^y c( b r) + 

+ s p a ^ ) n c ' r) +q3°' t) n°^ ) +qj“W‘p^ r) +q3f ) q3“«) 

(2? 3 n 6 w+2? 3 V 6 «)]. (104) 


Hence, the trace and the antisymmetric parts are com¬ 
pletely determined by the lower-order expressions of all 
parts (four equations for four components of P 1 -). What 
is missing is an equation for the five components of 
the trace-free symmetric term n* (rl) . Recall that the 
gravitoelectric system is only closed after imposing the 
MR, which then couples the trace-free symmetric time- 
evolution to the one of the trace and encodes the spatial 
dependence in a Poisson equation. 


B. Reconstruction of GR solutions 

To illustrate the scheme for the GR case, we will dis¬ 
cuss here how to reconstruct the full nth-order solution 
from the recursive equations {(94), (95)}. 


1. Trace part 

The trace part is the main part that is given by the 
hierarchy. In the absence of the trace-free symmetric 
term n* ( " ) , there is no antisymmetric term emerging and 
we are left with a recursion relation for the trace: 

p iP (n) = _l P M V2 p(r) _ Ip(‘)pMp 3 pM . (105) 

3 9 


2. Antisymmetric part 

It may appear counterintuitive that a nonvanishing 
antisymmetric part arises (starting from second order), 
given our assumption of irrotationality due to the given 
foliation of spacetime. However, this fact is known from 
the Newtonian Lagrangian perturbation theory, where 
antisymmetric parts arise, starting at second order, in 
Lagrangian space, while no vorticity is created in Eu- 
lerian space [15]. Our comoving setting corresponds to 
the Lagrangian picture of fluid motion, and the antisym¬ 
metric terms at order n satisfy and follow from the ir¬ 
rotationality condition (94), given all subleading terms 
p = 1... n — 1. However, we need to reconstruct a part 
of the trace-free symmetric term to recover all the New¬ 
tonian modes that have antisymmetric components, a 
problem to which we turn now. 
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3. Trace-free symmetric part 


the solution for the trace found from (87) reads: 


As our scheme does not separately provide a rela¬ 
tion that determines the five coefficients of the trace-free 
symmetric part (these equations are part of the gravit- 
omagnetic scheme), we have to reconstruct the relevant 
part that complies with the Newtonian solutions. To 
achieve this it suffices to realize that the one-form fields 
V a ( y X k , t) become integrable in the MR, d P l (X k , t), and 
so also the trace-free symmetric part. Hence, in the MR, 
the trace-free symmetric part of d P l (X k ,t) inherits the 
time-evolution from the trace. With this in mind, and 
due to the superposition property of our solution scheme, 
we are entitled to split the general trace-free symmetric 
coefficients into a part that reproduces the trace-free 
symmetric part of P a (X k , t ) in the MR, denoted by E , 
and another part H Rij- This is possible at any order: 

= Er , W £ n,' n ' , »(x t ) + H n ^ m \x k ,t). 

I m 

(106) 

The temporal coefficients £ <n,!> are the same for the trace 
and the trace-free symmetric gravitoelectric parts. For 
the full GR solution there is in addition a contribution, 
denoted by H Rij, which is related to gravitational waves. 
We will investigate this part in the followup article. For 
the time being we note that the superposition prop¬ 
erty discussed above assures that the resulting individual 
terms in the decomposition (106) are correct, if we use 
only this gravitoelectric part of the trace-free symmetric 
tensor in the hierarchy. Thus, even though the scheme 
does not determine all the components of P'j" 1 without 
solving the gravitomagnetic equations, it is consistent for 
the terms it delivers. Moreover, by inspection of corre¬ 
sponding perturbation and solution schemes that we de¬ 
rived for the gravitomagnetic part [3], we can conclude 
that the so-reconstructed solutions provide the leading- 
order modes of the relativistic solutions at any order. Of 
course, inserting the reconstructed solution into the full 
set of Einstein equations will result in constraints on ini¬ 
tial data in addition to the standard constraints. As an 
example we will discuss the constraints in the first-order 
scheme given below. 


C. Example 1: recovering parts of the general 
first-order solution 

To illustrate the hierarchy we begin with the first-order 
equations of the scheme {(94), (95)}, i.e. in coefficient 
form (87). With the split in space and time coefficients, 
the latter are the well-known solutions of the equation 
(equivalent to the equation in the Newtonian scheme [7, 
9, 10]): 

£ + 2Hi-^H i 2 n im a- 3 Z = Wa- 3 . (107) 

For an EdS universe the modes are proportional to a, 
a~ 3 / 2 and a 0 . Together with the initial conditions (51), 


p( i) 


3 

5 


(Uti + -Wt- 2 )a 


(Uti - Wti 2 ) a~i - *-WU 2 
(108) 


The antisymmetric part vanishes in view of (87), = 

0. We then need to reconstruct the trace-free symmetric 
part along the lines described in IV B 3 to complete the 
solution: 


E n% = ^ ( “imt i + 


3a 


u 


E T jtl 


5 a 3 / 2 


( E u%t i - E W%t i 2 ) - 


>E 


wit- 2 


(109) 


The notation tl stands for the traceless part. The initial 
fields have been split accordingly: 


Uij =: E U ZJ + H Uij ; Wij =: E W i5 + H W tj , (110) 


i.e. a part initializing the gravitoelectric, and the gravit¬ 
omagnetic parts, respectively. 

We remark that in Newtonian theory the tidal tensor 
is written in terms of the gravitational potential $: 

- £ij = - ^V 2 $ , (111) 

where a comma denotes derivative with respect to Eu- 
lerian inertial coordinates. If we consider the first-order 
solution (here restricted to the growing mode solution for 
notational ease), 

E F i j m = \w\t- 2 (a - 1) , (112) 

the first-order gravitoelectric part of the spatially pro¬ 
jected Weyl tensor assumes the form (note the conven¬ 
tional sign difference of this geometrical definition with 
the Newtonian (active) definition of £ij): 


fv 1 ’ = - fry 11 - 2#ny i) 

= -fi 2 Cd + 2Ha) 

= - § ( ' 2 “ (<V ~ l wsi 0 

= - h ( w '> ^ i ws ‘>) • < 113 > 

We find 

Eij w = - (w^ - ^ W6i jJ . (114) 


The trace W does not derive from a potential due to 
nonintegrability of the field. After executing the MR, 
we obtain (up to the conventional sign difference), the 
Newtonian tidal tensor (111). 





13 


Summarizing: given the formal analogy of the solution 
schemes discussed in Section HID, the above solution 
solves the gravitoelectric part of the corresponding rel¬ 
ativistic equations (87). The trace-free symmetric part 
(109), however, is only a part of the solution in the rela¬ 
tivistic case. Equation (102) in Ref. [LI] states that the 
first-order equation for the relativistic trace-free sym¬ 
metric part reads: 

= -a -2 , (H5) 


where is the trace-free part of the initial Ricci tensor. 
Plugging Eqs. (108) and (109) into Eq. (115), we can 
check whether our relativistic generalization satisfies the 
full equation. Three modes appear in the equation: a -2 , 
a -1 and cT^ 2 . The equation has to be satisfied at any 
time, and thus each mode must lead to cancellation of 
the coefficients. This leads to the following constraints 
{Hi := 2/3ii): 


% = - E W. 


Ejj\k _ Ejj^ 


(■ ij) \ k 


= a XJ 


k\ij i 


(ij) 

_ E 


(ij) ’ 


(ij) |k 


= E W 


k\ij 


(116) 


The first equation corresponds to the definition of the 
trace-free part of the initial Ricci tensor, Eq. (53) in the 
EdS case studied here. In view of the constraints U y 7 i = 0 
and W[jj] = 0 (c/. (50)), the other two conditions are 
equivalent to: 


E rr\k 


ur,,= E u 

ij | k 


k\ij 


E w ' [k ,, = E w k 


ij | k 


k\ij 


(117) 


What we call gravitoelectric part in the decomposition 
of initial conditions (110) is therefore determined to be 
the one that solves (116). The part contributing to 
the propagating gravitomagnetic part is then its com¬ 
plement. This labelling is not completely unambiguous, 
because in this scheme, the gravitomagnetic part com¬ 
puted from the gravitoelectric part is not null, see below 
for the first-order scheme. (Nevertheless, as we will show 
in the followup paper, it generates a null dynamical Ricci 
curvature tensor.) 

To check how constraining these relations are, beyond 
the constraints that we already have, we consider the 
first and second time-derivatives of the momentum con¬ 
straints and evaluate them at initial time in order to 
obtain constraints on the initial fields. Taking the sec¬ 
ond spatial derivative of these equations and contracting 
them with respect to one index, we get for Uij: 


U k j\k = U% 


Tjk _ Tjk 
^ j\ik ^ k\ij 


U' 


k\i __ jjk\i 


j\ik 


k\ij 


(118) 

The latter identity is solved by the gravitoelectric and 
the gravitomagnetic parts independently. For the gravi¬ 
toelectric part, we have: 


Erjk\i _ Ejj^i 

U j\ik ~~ U k\ij ! 


(119) 


which is equal to the once contracted spatial derivative of 
the above constraint (117). We conclude that Eq. (117) 
and the momentum constraints are compatible with but 
not equivalent to our constraints. They have to be solved 
independently in order for the solution to be compati¬ 
ble with both the evolution equation and the momentum 
constraints. 

What they do constrain are derivatives of the gravit¬ 
omagnetic part. To derive these constraints, let us first 
note that the first-order expression for the magnetic part 
can be found from Eq. (107) in Ref. [LI]: 

E Hij = a(ty\ E il jWs . (120) 


The solution for E Hij, cf. (109), shows that spatial 
derivatives of the first-order magnetic part can be traced 
back to spatial derivatives of Uij and W t j. Together with 
the first-order momentum constraints, P''i\ 3 = P l j\ii and 
imposing the constraints (116), we get: 


^uri ^si Ett n 


(i U j)s\lr 


-uri *1 E 


w. 


j)s\lr 


= 0 . ( 121 ) 


Thus, via (109), this leads to 

e uri E H ij , r = 0 , (122) 

i.e., the curl of E Hij vanishes. 

For its divergence the constraints (116) are not nec¬ 
essary. Taking the divergence of (120), and using the 
momentum constraints in the form E tS l l i is = 2/3P| is , we 
can show: 


E H ij \ i = 0 . (123) 

By combining (122) and (123), we conclude that 

A 0 E Hij = 0 . (124) 

Thus, the gravitomagnetic part that is generated by the 
gravitoelectric part is a harmonic tensor field at first or¬ 
der. This harmonic field can be constrained in the initial 
conditions (removed) by topological conditions on the 
perturbations. In an upcoming article, we will discuss 
such conditions. 


D. Example 2: constructing second—order solutions 
for ‘slaved initial data’ 

Let us now write out the system {(94), (95)} explicitly 
for n = 2. We simplify the first-order source by im¬ 
posing the so-called ‘slaving condition’ C7® - = (as 

explained in Refs. [9, 10], for second order in Ref. [11], 
and used to present the third- [12] and fourth-order [37] 
Newtonian solutions). This is not necessary but increases 
readability. The sum of (108) and (109) becomes: 

PV = lw ij t i 2 {a-l) . 


(125) 
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At second order (89) is simply 


2>i PT = - (v 2 Pf^ Pf> + (v 2 Pf ') P^ (1) , (126) 
and we have the system: 

f ? » + 2 £j<.. + 3|r = ^(«-‘-«-») i(127) 
[ c (2) = - ww, 

3 

with the source g <2> (f) = -ti 2 (a -1 — a” 3 ). 

To systematically determine the temporal coefficients 
of the hierarchy, it is useful to write the operator Pi in 
terms of a. We find: 


9 ( a ) — ^im^i X 

(G + “ 2c ) p " {a) + 1 + 2ac ) p,(a) “ i p(a) ) ’ 

(128) 

where c = For an EdS background, c = 0, the 

homogeneous solution is 

D (a) = aC\ + a~ 3/2 C 2 ; (129) 


Green’s function of Eq. (101) is 


G (a, s) 


2 s (a 5 / 2 - s 5 / 2 ) 
5 fl im Pi 2 a 3 / 2 


(!30) 


Now, it is a matter of a simple integration and Eq. (100) 
gives the second-order trace solution: 


p<2) = i C «2) a+ 2 C ,(2) a -3/2 + 9 ii 4 ^ + _ ( 131 ) 


To find the spatial coefficients of the solution, we use 
the initial values for the coframe and its time-derivative. 
They have been chosen to vanish for all orders higher 
than one in the hierarchy of solutions of Eqs. (89). There¬ 
fore, we find the system: 

ak 

P (2) (fi) = X C + 2 C + — C (2> fi 4 = 0 ; 

28 

P^\h) = |- X C - j 2 C + AC( 2 V = 0 , (132) 


which fixes all constants to be oc C (2) . Thus, the second- 
order trace solution reads: 


p(2) = £» (' WW _ W)W{) ; 


£ 


( 2 ) _ 
+ — 




-3/2 


. (133) 


After executing the MR this coincides with the second- 
order Newtonian solution of Ref. [15]. 

The antisymmetric equation (94) still delivers fp a j (2) = 
0. This is due to the restriction to ‘slaved initial con¬ 
ditions’, otherwise we would have a nonvanishing part 


here. Thus, we only need the trace-free symmetric part 
to complete the solution. The gravitoelectric part can be 
written as 


£ n< 2) = £V 2, 4 2) , (134) 

where the trace of is given by (\V 2 — t 2 . 

The rest of its components can be determined from 
the generalization S^fj where S ^ is the so¬ 

lution to the Newtonian Poisson equation Ao6>^ 2 * = 
i^W k \k) 2 — t 2 , and where A 0 denotes the 

Laplacian in local (Lagrangian) coordinates (see Ref. 
[15]). To avoid passing by the generalization of the New¬ 
tonian result, one can of course also insert into 

(134) and solve the remaining relativistic equations of the 
gravitomagnetic part to find the off-trace components of 
o b . The gravitomagnetic contribution encoded in <Sb 
has recently been studied to second order in Ref. [39]. 

The explicit derivation of the inhomogeneous second- 
order term in this subsection illustrates that, using (100) 
and (130), the calculation of the temporal evolution of 
the general relativistic trace part is straightforward and 
only involves the calculation of integrals. This can also 
be easily extended to perturbations of a ACDM universe 
model by noting that (129) becomes: 


D{a) = a (2) Pi Q,l,^;-c/) Ci + \J-^ + cC 2 , 

(135) 

with the Gauss hypergeometric function ^Pi. Greens’ 
function reads in this case: 


G(s,a) = |- S 


5 


D+ (a, c) - D+ (s, c) 


where D + (a, c) is the first term in (135). 


(1 + ca 3 ) s 3 
(1 + cs 3 ) a 3 

(136) 


V. SUMMARY AND CONCLUDING REMARKS 

We have investigated gravitoelectric perturbation and 
solution schemes at any order in relativistic Lagrangian 
perturbation theory. These schemes cover the full hier¬ 
archy of the Newtonian Lagrangian perturbation theory 
if restricted to integrable Cartan coframe fields. 

Despite the fact that the solution scheme presented 
in this work gives on its own not all parts of the rela¬ 
tivistic perturbation solutions, it delivers an important 
part relevant to the formation of large-scale structure. 
As is well known (see, e.g., discussions in Refs. [31] and 
[3]), the fastest growing scalar modes of the GR solu¬ 
tions correspond to the Newtonian modes, shown up to 
second order and, by inspection of the schemes we investi¬ 
gated, we showed this to hold for the gravitoelectric part 
also beyond second order. As we recover all the New¬ 
tonian terms with their correct temporal evolution and 
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their constrained spatial coefficients, we also know that 
our solution contains all terms that become important 
in the Late Universe. The presented scheme is explicit 
enough to derive solutions at any desired order by al¬ 
gebraic codes along the lines of the reconstruction rules 
that we exemplified up to the second order. 

We demonstrated the close formal correspondence of 
the gravitoelectric Lagrange-Einstein system to the New¬ 
tonian theory furnishing construction rules that also al¬ 
low to find other, nonperturbative relativistic solutions 
from Newtonian ones. 

The role of gravitational waves, corresponding to the 
missing part in our scheme, has to be further explored. 
The missing part, which we denoted by H Hij in the co¬ 
efficients of the trace-free symmetric parts of the pertur¬ 
bations, corresponds at first order to ‘free gravitational 
waves’, i.e. that part of gravitational radiation that does 
not scatter at the sources. This changes at higher or¬ 
ders, since this part will couple to the sources starting 
at second order. We will investigate in detail the gen¬ 
eral first-order scheme including gravitational waves in 
the next article of this series, where we also identify the 
transformations and restrictions that have to be imposed 
to obtain the known solutions of the standard perturba¬ 
tion theory, where perturbations are embedded into the 
background spacetime. 
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Appendix A: Equivalence of integrability of the 
coframes and the flatness of space 

The standard choice of orthonormal coframes fj a in the 
Cartan formalism implies for the spatial metric coeffi¬ 
cients gij = S a bfj a ifj b j, with rj“ (tf) ^ <5“ at initial time, in 
order to have an initially nontrivial metric. 

The alternative choice of adapted coframes r/ a , used 
in this article, represents the metric coefficients as g t j = 
Gabrf'i'n 0 'j, where we are entitled to require 77 “ (tf) = 5 a i 
at initial time, encoding the initial metric into the coeffi¬ 
cients Gab, he. = G a bS a iS b j. This makes the compari¬ 
son with the Newtonian choice of Lagrangian coordinates 
to coincide with the Eulerian ones at some initial time 
more direct. 

As we discussed in (36), the basic assumption is that 
both coframe types describe the same metric form, i.e. 


g = 8 cd f] c 0f] d = Gab^l ° ® T? h , from which we infer: 

Gab = S cd rj c a fi d b , (Al) 

where we denote with f) c a = efff\ the coefficients of the 
projection of rf onto the basis r\ a . 

The MR applied to either of these coframes requires 
them to be exact forms, f] a = d/° or r/ a = d f a . They 
then define some global Eulerian coordinates, x a and x a , 
respectively. In the MR, Eq. (Al) is equivalent to 

dfC Qfd 

(A2) 

We infer from (A2) that the coefficients G a b just depend 
on initial vector displacements after executing the MR. 
They are related to the initial deformation gradient in 
the orthonormal description, as can be seen by looking 
at the metric equivalence relation in an exact Lagrangian 
basis, g = g., j dX l 0 dX J = g^dX 1 0 dX J : 

g (X k ,U) = 5 ab f a ^ i {X k ,t i )f b ^ j {X k ,t l )dX i 0 dU 

= 6 ab f a li (X k ,t i )f b lj (X k ,t i )dX i 0 dXi 

= G a6 (A fc )<5“ ( 5 h J dA i ® dX* , (A3) 

where a slash denotes derivative with respect to the coor¬ 
dinates X ', as in the main text, and it is explicitly noted 
otherwise. From (A2) we conclude that 

df c df d 

g = G ab r] a 0r] b = 8 cd -^-^dx a 0dx b = 5 cd dx c 0dx d , 

(A4) 

which is the Euclidean metric. 

Summarizing: execution of the MR leads, in either of 
the chosen coframes, to a metric that is equivalent to 
the Euclidean metric. The coefficients G ab can then be 
expressed in terms of initial vector displacements, cf. Eq. 
(A3). 

Appendix B: Erratum and Remarks 

We correct a mistake in the paper [20], add a clarifi¬ 
cation to the paper [LI] concerning the gravitomagnetic 
part of the spatially projected Weyl tensor. 


1. Newtonian perturbation scheme, Ref. [20] 


The perturbative system of equations derived from La¬ 
grangian Newtonian theory leads to {(65),(66)}: 

% d P i A (dX :l + d P J ) = 8 ija - 2 dlP A dX> ; (Bl) 


tijk 


((D - 4nG q h ) dP k ) A dX j A dX k 


+ ((2 D - 4ttGq h ) dP l ) A dP j A dX k 


D - -Tp-QH ) dP* ) A dP j A dP fc 


— Eijk~ 


47 tG . 


Sgia-MX 1 A dA J A dX k , 


(B2) 


3 
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where the operator D := d 2 /dt 2 +2HcL/dt. The numerical 
coefficients of this system differ from the result presented 
in Ref. [20]: Eqs. (28a-d). 

2. Gravitomagnetic part of the spatially projected 
Weyl tensor, Ref. [LI] 

In Ref. [LI] the symmetrization of the spatial parts of 
the Weyl tensor has not always been written explicitly, 
which may lead to confusion. The idea of a not manifestly 
symmetric writing is best seen in the equations for the 
gravitoelectric part: its definition in Eq. (32) already 
assumes that the field equations hold, which are then 
recovered by the conditions (33). While having advan¬ 
tages, this representation is implicit. The same applies 
for the gravitomagnetic part (see Ref. [LI]: Eqs. (70,73), 
where the momentum constraints have been inserted and 
recovered through an explicit symmetry condition). If 
we wish to consider the original geometrical definition of 
these tensors without inserting the field equations, then 
we have to write, e.g. for the gravitomagnetic part (here 
written for the adapted coframes used in this article): 

Hij I \l 

= ~ 27 £ "“ 9 "< iG “ l ((’ , “>>’ i ‘*)||, + (*’b>)||,) 

= - 7fJ enkl 9n(iGab (v a m v\ + 

+ r l a k\\iV b J )+V a kV b M i) , (B3) 

where g n i = G a bV a n V b i • Employing the irrotationality 
condition (28) and taking its covariant derivative, 

GabV a {i r l b j)\\k + G a bV a \jV b i]\\k = 0 > (B4) 

we obtain an expression for H tJ that manifestly vanishes 
with the passage to the Newtonian theory through the 
MR: the covariant derivatives of the (now integrable) 
coframes vanish: 

Hij = --j£ nkl gn(iGab (jl°j)tf= 

1 nkl r i /~i fd re ra ( rb rb i m re \ _ r\ 

'-'de'-'abj \nJ \(il \j) \yJ \kl ~ J \m n ,c J \kl J ~ U ) 

(B5) 


which is the result claimed in Ref. [LI]: Eq. (84), and 
which is expected in a concise execution of the Newtonian 
limit [21]. 


3. Remark on the Newton—GR “dictionary” 

The reader may have noticed that the Newtonian ten¬ 
sors which correspond to the GR tensors in the MR dis¬ 
play mixed indices, e.g. in Subsection IID 3 for the MR 
of the gravitoelectric part E l 3 . We explain here why this 
is the case. 

Consider, e.g., the Newtonian field strength gradient, 
g a b , where a comma denotes a derivative with respect 
to Eulerian coordinates. We have on purpose denoted 
the vector index by a counter index here, counting the 
number of vector components, and we have also used 
the counter index for the Eulerian derivative, since both 
become noncoordinate indices by executing the inverse 
MR, i.e., by passing from the Newtonian deformation 
gradient with respect to Lagrangian coordinates, /“ fe , 
to the nonintegrable coframe coefficients, r\ a k . This ex¬ 
presses the fact that the vector embedding space disap¬ 
pears through this operation; in the integrable case, the 
components /“ can be considered as coordinate functions 
with the coordinates x a “W Now, the transformation of 
the field strength gradient to Lagrangian coordinates in¬ 
volves the inverse of the transformation /“, which we 
denote by X 1 = h l (x b ,t): g a b = f a ^h l b , where we have 

also inserted the definition g a = f a . Denoting the co¬ 
efficients of the inverse of the coframes by e b , we see 
that the inverse MR constructs the relativistic analog to 
g a b , which is T a b := ri^e^. Projecting this latter field 
onto the exact basis with the help of the coframes and 
their inverse (the frames), we obtain a field with mixed 
indices, T a b e^r] b j = T' 3 . The tensor coefficients are then 
obtained by lowering an index with the spatial metric, 
•Rfei = guFj = G ab r] a k i) b 3 . For mixed indices, we obtain 
(compare Eq. (46)): 

r j = e a% = ^e a6c e iH rj“ ■ (B6) 
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